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Abs tract ~Th.e  theoretical  analysis  of  an  enclosed  air-bridge  is  presented.  This 
includes  a  derivation  of  the  Green’s  function  in  the  x  and  z  directions  which  is  used 
to  find  the  current  distribution  on  the  conducting  strips  of  the  air-bridge.  Boundary 
conditions  at  the  interfaces  are  applied  and  the  numerical  technique  Method  of 
Moments  is  used  to  solve  the  integral  equation  for  the  unknown  current.  Upon 
derivation  of  the  current  distribution  on  the  conductors,  an  ideal  transmission  line 
model  is  applied  to  obtain  the  scattering  parameters  of  the  structure. 
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1  Introduction 


Millimeter  wave  technology  concerns  itself  with  that  portion  of  the  electromag¬ 
netic  spectrum  between  0.3  Ghz  and  300  Ghz,  corresponding  to  wavelengths  of 
1000  mm  to  1mm.  Effective  quasi-static  techniques  have  been  developed  for  the 
lower  frequencies  (0.3  Ghz  to  3  Ghz)  but  for  the  higher  frequency  part  of  the 
spectrum,  a  full-wave  analysis  must  be  employed. 

Millimeter  and  microwave  systems  may  be  overshadowed  by  infarred  and  optical 
systems  but  limitations  to  the  latter,  in  particular  their  disadvantages  in  fog,  dust, 
rain,  and  nighttime  viewing  support  further  development  of  the  former.  As  with 
many  technologies,  the  number  of  applications  will  increase  with  the  passage  of 
time. 

The  typical  millimeter  microwave  integrated  circuit  contains  associated  active 
and  passive  elements  interconnected  by  transmission  lines.  In  integrating  these 
components  together,  various  discontinuities  arise  where  evanescent  fields  and  sur¬ 
face  waves  play  an  important  role  in  their  operation.  Work  has  been  done  to  model 
these  discontinuities  with  lumped  elements  but  the  numerical  techniques  used  to 
derive  equivalent  circuits  are  either  frequency  bound  or  dependent. 

Here  we  present  a  full-wave  analysis,  which  in  not  frequency  bound  or  depen¬ 
dent,  to  analyze  an  air-bridge  structure.  The  resultant  expressions  and  methods 
used  are  general  enough  to  be  applied  to  an  array  of  three  dimensional  problems. 
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a.  On*  dtmanwonal  view  nignugnting  curort  airaetton* 
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c.  Thr**  dhnanatonal  vM*t  NgNIghllng  raiatfr*  gaomatiy 

Figure  1:  Air-bridge  in  enclosed  microetrip. 


2  Evaluation  of  the  unknown  current 


In  our  problem  the  current  distribution  on  the  structure  of  interest  must  be 
determined  accurately.  Then  by  the  use  of  an  ideal  transmission  line  model,  the 
scattering  parameters  can  be  evaluated.  In  obtaining  the  current  distribution, 
Pocktnigton’s  integral  equation  is  solved  numerically. 

The  formulation  of  the  Pocktnigton’s  integral  equation  and  the  solution  for  the 
given  structure  are  presented  in  the  following  work.  Our  structure  under  study  is 
shown  in  Fig.l. 
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2.1  Formulation  of  the  Integral  Equation 


Through  the  manipulation  of  Maxwell's  equations  1 

VxT=  —jjjjiR  (l) 

V  x  77  =  jufiE  -t-  J  (2) 


V  •  tE  —  p 

V  ■  nH  =  o 


along  with  the  representation  of  the  magnetic  vector  potential  A 


77  = 


Vx  A 
(* 


(3) 

(4) 

(5) 


one  arrives  at  an  expression  relating  current  and  the  magnetic  vector  potential 


V2A  +  k‘A  =  -tiJ. 


(6) 


When  the  current  7  is  represented  by  a  dirac  delta  function  in  equation  (6),  the 
Green's  function  becomes  a  solution  as  shown  by  the  equation 

V3£ +  **£:=  (7) 


To  obtain  a  unique  solution  that  applies  to  the  specific  geometry  as  shown  in 
Fig.l ,  one  must  apply  the  characteristic  boundary  conditions  of  the  structure.  Note 
that  we  have  introduced  the  dyadic  form  as  we  need  to  be  able  to  describe  fields 
which  are  produced  by  a  current  of  arbitrary  orientation. 

'throughout  ihtt  report  «a  eJ*‘  tiam  convration  it  tttumtd  to 6  mpprwid 
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For  the  case  of  a  single  x-directed  current,  the  unit  dyadic  I  takes  on  the  form; 

7  =  xx.  (8) 

The  vector  representation  of  this  current  is 

J  =  6(r  -  r')x.  (9) 

Equation  (9)  represents  a  dipole  directed  in  the  x-direction  and  parallel  to  the 
interface  between  regions  II  and  III.  It  has  been  shown  by  Sommerfeld  that  the 
magnetic  vector  potential  of  this  structure  needs  to  have  two  components  so  that 
the  appropriate  boundary  conditions  are  satisfied.  This  dictates  that  A  must  have 
one  component  parallel  to  the  current  source  and  another  parallel  to  the  interface; 

7t*  *Aix  +  A\z.  (10) 

The  integral  equation  which  relates  the  magnetic  vector  potential  to  the  current 
of  interest  is  written  as 

X  =  nj  f  Jy7?  7 dv  (11) 

where  the  dyadic  Green’s  function  7T  is  uniquely  defined  by  the  structure  under 

investigation  and  takes  on  the  general  form 

G'^xy  G\,xz 
G\xyi  G'yvyy  G\tyz  • 

Q^zy  Gi„zz 

]i  denote*  for  which  region  (I,  II,  or  III)  the  equation  appliae  «•  defined  in  Fig-2 
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In  the  case  of  a  single  directional  current  in  the  x-direction,  the  dyadic  equation 
has  only  two  components, 

(f  =  cr„xx  +  (12) 

From  equations  (l)-(5)  our  electric  field  is  related  to  the  magnetic  vector  po¬ 
tential  by 

V  xH 


therefore 


r  =  =  — V  x  (V  x  A)  =  —{k2A  +  W  •  A) 

jut  jutp  jutfi 


.  f.  1  6  ( 6At 

Ex  [A‘  +  (**)2  Sx\Sx  +  Sz)\ 

_  .  1  6  (A,  *A,\1 

Ey  JU'[(*<)2*y  \6x  +  6z  jj 


E,  =  —ju 

and  through  Maxwell’s  equations, 


,  ,  1  6  (6A,  6A.  V 

*+  (M2**  \  6x  +  J 


Hr  = 


16A, 
H  Sy 


L/sa1_sal\ 

fi  \  OZ  OX  J 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


2.2  Derivation  of  Green’s  function 


For  our  structure  pictured  in  Fig.l  (an  air  bridge  on  an  enclosed  thin  microstrip), 
we  require  only  the  derivation  of  the  Green’s  function  for  an  x-  and  z-directed  cur¬ 
rent.  For  the  variation  of  current  in  the  y-direction,  we  have  assumed  a  Maxwellian 
distribution.  Our  first  step  will  be  to  derive  the  Green’s  function  considering  the 
x  and  z  components  of  the  current  separately. 
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Figure  2:  x-directed  current  above  dielectric  in  enclosed  microetrip  structure. 

2.2.1  Derivation  of  the  Green’s  function  for  a  x-directed  current 

With  equations  (14)-(19)  the  Green’s  function  can  be  related  to  the  electric  and 
magnetic  fields  which  have  conceivably  defined  values  dictated  by  the  structures 
electrical  characteristics.  We  begin  by  stating  that  the  tangential  electric  fields  are 
zero  everywhere  on  the  surface  of  the  walls  of  the  structure; 


E*yt  =  0  at  x  =  0,  a 

(20) 

E'ta  =  0  at  y  -  0,  b 

(21) 

Elx<y  =  0  at  *  =  0 

(22) 

Ex,l  -  0  at  z  =  c. 

(23) 

Employing  separation  of  variables  to  the  expressions  and  the  established  bound¬ 
ary  conditions  (20)-(23),  the  following  general  forms  of  the  green's  functions  can 
be  derived  for  each  region  of  interest 
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G*x  =  £  5Z  ^cos  f - \  sin  (— 3in{k[z)  (24) 

n=Om=l  \  a  /  \  0  / 

G'„  =  ££b^»(=E)«.(22)«.(*J,)  (25, 

n=lms|  \  a  /  \  O  / 

G'i  =  52YlC09  (““)  s*n  (Ansm(kllz)  +  C!Icos{klt!z))  (26) 

G"  =  £  £>s  (^)  sm  (^p)  (*"«»(*",)  +  D“cos{k['z))  (27) 

nsOmsl  \  a  /  \  0  / 

<?"'  =  ££A"W^)*(^Wtf"<.-e))  (28) 

nsOmsl  \  a  /  \  0  / 

G'”  =  t  t  (^  )  sin  (^p)  3m(*"'(,  -  C).  (29) 


In  equations  (24)-(29)  the  eigenvalues  ib„  kv,  and  kt  satisfy  the  following  rela¬ 
tions: 

(*•)’  =  (*;>' + (*;)’ + (*;>a 
where  it,  = 

»nd  k,  =  (Sp) 

k 2  *  u »V«.  (30) 

To  determine  the  unknown  coefficients  AUI,in  and  BIuin ,  one  should  apply 
boundary  conditions  at  the  two  interfaces  between  Regions  I,  II,  and  III.  Since  at 
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the  boundary  between  Regions  I  and  II  there  exists  no  magnetic  charge,  the  normal 
magnetic  field  must  be  continuous  across  the  boundary.  Therefore  if/  =  if/7,  so 


Alsin{kiz')  =  An  sinik*1  z  )  +  Cn sinik1/  z)  (31) 

Since  there  is  no  electric  current  in  the  y  direction,  the  tangential  magnetic 
field  is  continuous.  By  substituting  this  relation  into  (18),  one  obtains  H[  —  . 

therefore 

Blcos(k[  z)  =  Bnsin{  ki'z)  +  Dncos{] IfcJV)  (32) 

Also,  since  there  is  no  magnetic  current  in  the  x  direction,  the  tangential  electric 
field  is  continuous.  Therefore  Ely  =  £7/; 

A^in^z')  —  Blsin{klxz)  — 

B11  cos{k[I z)  -  Dtlsin(kitx)  +  Allsin{k{lz')  +  C"«n(*"z')  (33) 

By  applying  similar  boundary  conditions  on  the  interface  between  Region  II 


and  III,  the  following  equations  result; 

Ansin(kllH )  +  Cncos{kl/H)  =  AlIIsin(kin(H  -  c))  (34) 

Bu ain(k[l H)  +  DIlcoa( k'/H)  =  Butco3(k[ll{H  -  c))  (35) 

k[l (A11  cos{k[r H)  -  Cl!3in(kilH ))  =  A1"  k[u  cos{k[n  (H  -  c ))  (36) 

AIl3in(k[!H)  +  C"co3(klglH)  (37) 

( BlIco3{k[lH )  -  DlI3in(k[lH)) 

V  a  / 


=  —  (V"»m(t"'(ff  -  c))  +  -  c)))  . 

«///  \  \~T'  * 
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Since  there  does  exist  an  electric  current  between  Regions  I  and  II,  the  mag¬ 
netic  field  between  these  two  boundaries  is  discontinuous.  This  discontinunity 
cam  be  considered  by  integrating  the  inhomogenous  helmholtz  expression  over  the 
boundary  and  using  orthogonality.  This  results  in 

—  (kIzArcos(klzz)  —  kIzIAllcos(kIzlz)  +  k[lCn  sin(klJ  z')) 


-— (¥)•*(=£)  <38> 


where  p  — 


( 


\ 


4  when  n  ^  0 
2  when  n  =  0 


Using  these  eight  relations  derived  from  the  boundary  conditions,  the  eight 
unknown  constants,  A1,  Bl,  A11,  Bn ,  Cllt ,  D 11 ,  A111,  Bin  are  found.  The 
resulting  Green’s  functions  for  a  x-directed  current  above  the  dielectric  upon  sim¬ 
plification  are 


GL  = 


EE 

nsOmil 


abk,[cot(ktz')cos(kMz')  -  stn(ifcfz')j 


CO* 


39) 


'/  _ 


00  00 

=  EE 


/K1  -  tm’(t.M))  (?)  (<"'  -  l)tan(k',"(H  -  c)) 


ob[klirtan(ktH)  -  kMtan(k{n(H  -  c))][coi(ktz')co3{kzz')  -  sm(k 

1 


) 


cos 


(=?)*(=?)? 


k‘"tan{k{"{H  -  c))  -  el"k,tan(k,H) 
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/ nirx\  .  /mny\  . 

sin  [ - 1  stn  ( — —  J  cos(kxz) 


(40) 


^  ~  _ p  [kj11  -  kttan(kjn(H  -  c))tan(kzH) 

n=oi  «6fc.[arf(fc,*')cM(Ar,*')  -  5in(*,z')][*'"fan(fcfff)  -  kztan(k*Il{H  -  c))] 

MJ  (^r) sin  Or*)  «"<*."*)«»  (^p-) «» (:2p-)  + 


aMrr[c<?£(£,2,)cos(/bI.z')  —  5in(A:,z')]< 
/  nirx\  (  miry'\ 

c°s{—)3"'  (— j 


=  EE 


P  (1  -  tan\k.M))  (m)  (t"1  -  l)t<m(*'"(g  -  c)) 


toZZi  ab[k{ntan(ktH)  -  k,tan{k[n{H  -  c))][cot(fc,i:')c<w(fc**')  -  sin{ktz')] 


1 

f  nirx's 

1  Qi’n  1 

rmiry\ 

—  c))  —  elntan(kxH) 

1  3,71  ! 

l  b  ) 

( nirx'\  .  fmiry'\ 


~  V'  V' _ P  (1  —  tan2{ktH))cos(kzH ) _ 

ho  hi  ab[cot(kgz)cos(ksz')  -  sin{kxz')]cos{klxn{H  -  c)) 

1 

[klx"tan{k,H)  -  kxtan{k["(H  -  c))] 

cos  (52)  Jin  (2p)  -  c))a»  (2^)  sm  p 


”  ”  p(l  -  tariff))  (?)  (el"  -  l)tan(fc"'(//  -  c)) 

ab[kl"tan(ktH)  -  kxtan(klxll(H  -  c))][cot{kxz')cos{kxz‘)  -  sinf*.*')] 
sec(kl"(H  -  c )) 

[k'"tan{k,(H  -  c))  -  e'"*an(Jt'"ff)] 


(?) s,n  (?)  cw<*"'<z  -  c»«*  (?) sm  (?) 
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Figure  3:  Structure  when  conducting  strip  is  lowered  on  to  the  dielectric. 


If  the  conducting  strip  is  lowered  on  to  the  dielectric,  we  will  have  only  two  regions: 
air  representing  Region  I  and  dielectric  representing  Region  III. 

To  formulate  the  following  Green’s  fuctions  for  the  strip  on  the  dielectric  ,  one 
simply  lets  z  —  H  in  the  previous  expressions. 


P  tan(k{lI(H  -  c)) 


5,5  <‘Hcos(k,H))lk‘"Um(k,H)  -  -  c))j 

(45) 


(mxy\ 

COS{  a  )  1 

K  b  ) 

,  _  »  P  tan(k,H)  («)  (<"'  -  1  )tcn(kj"{H  -  0) 

C"  ”  “oir,  *b(coi(k,H))lk<"tan(k,H)  -  k,tan(k‘,"(H  -  c))l 


1 


[k[utan(k[Il(H  —  c))  —  tlTn  kttan(kxH)]C°3 


G','J  =  EE 


p  tan(ktH) 


ZZotZx  abcos{k["{H  -  c))[k["tan{k,H)  -  k,tan( k“‘(H  -  c ))] 

(=r)  -  (T)  *■«"<*  - fr1)  -  (=?>> 


»///  _ 


=  EE 


P(4“  -  l)t<"»(t,'"(g  -  c)t<m(t,tf)  (?) 


abco*(k,(H  -  c))[k<“tan(k,H)  -  k,tan(k‘"(H  -  c))] 

[It ["tan(k‘,‘<(H  -  c))  -  C‘"k,tm(k,H)\‘in  ("r)  "B  (^P)  ~ 

Cnirx'\  f  mir  y'\ 

—)”n{—) 


c)j 


(4S) 


2.2.2  Derivation  of  the  Green’s  function  for  a  i-directed  current. 

The  structure  used  to  consider  the  z-directed  current  is  pictured  in  Fig.4.  We 
will  simplify  our  problem  separating  it  into  two  parts;  a  primary  field  problem 
and  a  secondary  field  problem  as  shown  in  Fig.5.  For  a  z-directed  current  only 
one  component  of  the  magnetic  vector  potential  is  needed  to  satisfy  the  boundary 
conditions; 

A  =  zAz  (49) 
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Figure  4:  A  z-directed  current  above  a  dielectric  substrate. 


X 


z 


z 


Figure  5:  The  problem  ia  divided  into  two  parts,  a  primary  and  secondary  problem. 
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where 


a;  =  a'„ + a;. 

a,"  =  A"  +  A" 

A,'"  =  A"'. 

From  our  definition  of  Magnetic  vector  potentials, 

H  =  hl 


(50) 

(51) 

(52) 


(53) 


and  from  Maxwell's  equations, 

Vxl 


E  = 


jut 


-A- V  X  (V  X  A)  =  r-i~(fc2A  +  W  •  A) 

JUtfjl  ju/e^ 


therefore,  for  the  z- directed  current, 


(54) 


,_L(*>+£U 

june  \  oz 2  ) 

(55) 

(56) 

(57) 

ii 

o 

(58) 

»-H#) 

(59) 

(60) 

Now  we  must  apply  the  boundary  conditions  of  the  structure  to  obtain  the 
general  (primary)  solution.  The  first  readily  known  boundary  conditions  are  those 
on  the  walls  where  the  tangential  electric  fields  become  zero.  Therefore 


Ei  =  0  at  z  =  0  y  =  0, 6 


(61) 
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(62) 


£"  =0aty  =  0,b 

£'  =  0  at  z  =  0;  x  =  0,  a  (63) 

=*  0  at  z  =  0,a  (64) 

£f  =  0  at  x  =  0,a  y  =  0,6  (65) 

E[l  =  0  af  x  =  0,  a  y  =  0, 6  (66) 

In  the  primary  field  problem,  the  inhomogenous  differential  equation  takes  on 
the  form 

VM,  +  k20A,  =  -fi.  (67) 


The  solution  to  the  above  equation  when  the  listed  boundary  conditions  are  satis¬ 
fied  is  of  the  form 


I,  -  £  £  A‘>in  ( 

'nrx\  .  / 
—)  stn\ 

:=f) 

cos(fcfz) 

(68) 

nail  mat  ' 

oo  oo 

"  = 

nal  mmO 

{—)  Sm 

(miry' 
\  b  , 

(69) 

In  the  secondary  field  problem,  we  will  derive  a  solution  that  satisfies  the  ho¬ 
mogenous  differential  equation.  This  is  due  to  the  fact  that  we  do  not  have  a 
current  source  in  the  secondary  field  problem.  In  both  field  problems,  the  electric 
fields  must  satisfy  the  same  boundary  conditions  on  the  conducting  walls.  As  a 
result,  the  secondary  fields  are  of  the  same  dependence  with  respect  to  the  x  and 
y  coordinates,  therefore 
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(70) 


A"  = 


£  £(/"-*m(M)  +  fjlcos{kzz )) 

nsl  mxl  V  / 

.  (nirx'\  .  ( miry'\  , 

in  I  —  I  sin  I — - —  1  cos[ktz  ) 


sin 


(71) 


For  Region  III  we  have  the  same  standing  wave  solution  as  in  Region  1  except  for 
the  fact  that  the  conducting  wall  has  been  moved  by  (H+h)  along  the  z-axis  as 
reflected  in  our  z  dependence  below 


w  w 

Al"  =  EE; 


4  ft 


mjk,z  flit 


■■  (¥)■'•(=?) 


in  (^P)  sin  (^p)  cos(kin(z  -{H  +  h)). 


sin 


(72) 


In  the  case  of  a  delta  function  source,  Atp  will  give  a  G„p  component  in  the 
dyadic  Green’s  function  and  similarly,  a  A„  will  give  a  Gtxt.  These  two  components 
are  related  to  G„  by  the  relation 

Gtt  =  Gttp  +  Ggt t  (73) 

Since  we  have  a  current  source  between  Regions  I  and  II,  our  magnetic  held  is 
discontinuous  so  we  must  integrate  the  inhomogenous  helmholtz  equation  over  the 
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interface  which  results  in 

lim  /  (V2  +  k\)Atdz  =  ~n  lim  [  6{x  -  x)6(y  -  y')6{z  -  z')dz,  (7-4) 

—a  a**0  Ji  —a 

and  upon  simplification,  one  obtains 

]»n  (J;A*  =  -/**(*  -  x')6(y  -  y)  (75) 

f  All  CAl 

~ir  “ir  '«'■  ^(x  -  *  w*  - *>  (76> 

from  (25),  one  obtains 

J  (kisin(kiz')Ai  -  jk[l  A")  -  sin  «n  •  (77) 

One  more  equation  is  needed  to  solve  for  the  two  unknowns  A£,  and  A™ .  Utiliz¬ 
ing  the  fact  that  the  E  field  in  the  z-  direction  must  be  continuous  at  the  boundary 
between  regions  I  and  II;  one  obtains 


Alcos(ktz)  —  A11 . 


(78) 


These  equations  are  solved  for  the  unknown  coefficients  resulting  in 

<  =  (^7-)  Vy 

By  substituting  these  expressions  into  our  general  forms  (68)  and  (69),  the  zz- 
component  of  the  Green’s  function  takes  the  form; 

«•  - 


stn 


(81) 
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sin 


(^)-(T)'-'‘! . ' 


(82) 


In  order  to  determine  the  four  unknowns  of  the  secondary  held  problem,  bound¬ 
ary  conditions  on  the  air- dielectric  must  be  applied  to  the  structure  of  Fig. 5c. 

First,  the  z-component  of  the  electric  field  is  continuous  across  the  boundary 
between  regions  I  and  II;  £/  =  E[l .  Therefore 

fIcos{kzz)  —  fjlsin(ktz)  —  flJ cos(kzz  )  =  0  (83) 

Integrating  the  homogenous  helmholtz  equation  (there  is  no  electric  source  in 
the  secondary  problem)  across  the  boundary  and  using  orthogonality,  one  obtains 

///cos(fcJz')  -  fl1 sin(kxz' )  +  fIsin(kzz')  =  0.  (84) 


From  (83)  and  (84),  we  conclude  that 


/."  =  0 


(85) 


and 


fl  =  //'• 


(86) 


From  the  boundary  conditions  at  the  dielectric  interface,  Ez!  =  Ezu  one  obtains 


je~^h*Hcos(kaz  )  +  fj!sin(ktH)  =  -  n^ui~  ^ 

cr  r*r 


(87) 


and  from  H!tl  =  HlxI! , 


e~}h,Hco3(ktz')  +  f'tcos(kzH)  1  euiCos(k[^h) 

lgt  ~  „///'  liii 


k{" 


(88) 
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The  resultant  Green’s  function  for  this  problem  will  have  a  zz  component  only, 
which  is  given  by 


G[,  = 


tt  +/')>&«»(=£)»» 


n=l  msl 

sin 


( m*y" 
\  *>  . 


Gil  = 


/niri\  / miry\  /t  x 
in  I - J  sm  ( — —  j  cos(kzz) 


H  £  j~hTelk‘Z  (e'jk,,cos{kxz)  +  firsin{ktz)  +  f[l cos{kxz)) 

n=l  m=l  a°k* 

(nirx'\  .  (miry'\ 

—)3m  [—) 

-■CtM?) 


(89) 


(90) 


Gi"  = 


oo  oo 


4^ 


-  (=?)*(¥) 


V  V  i  M  //// 

sin  OrO 3in  (tit)  cos(k*n(z  ~(H+ h )))• 


(91) 


After  applying  all  the  necessary  boundary  conditions,  the  Green's  function  takes 
on  the  following  form 


r/  _  A  c^co^M*  -  ff))cos(fc'"h)  4-  -  H))sin(k“‘h) 

”  mIimi  k*elnsin(ktH)cos(klxl!h)  +  klxIlkxsin(klxIlh)cos{kxH) 


4 p  .  ( nrx  \  .  ( mxy  \ 
a  )3m\  b  ) 

sin  (22)  am  (22)  «*(*.*) 


(92) 
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G 


ti 

XX 


y.  Y  4"k.cos(k,(z  -  H))cos(kinh)  ±  k[u sm{k,(z  -  H))sin{k[11  h) 
^erm3zn(^/f)c03(A:^^)  +  kt"k23in{k["h)cos{ktH) 


4/z  (nirx'\  .  (miry'\  , 

(~  J  ,,n  (“  J CM(*-2  > 


/nirx\  ( miry\ 
Stn  {—)  3m  {—) 


(93) 


G{{1  =  ££ 


e^7/ kznl.n  cos(kzz ) 


SA  kje^Isin(ktH)cos(kIzIIh)  +  klzIl3in(k[uh)cos{ktH ) 
4/z  .  (nirx'\  .  frmry'\  < 

_.,ln  (— J  Jin  j  cos(*,z) 

«»  (“t9  sin  (^P)  c°s(ki"(z  -  (* + ^))) 


(94) 


2.2.3  Summary  of  Green’s  function  determination 

In  this  chapter  we  have  determined  the  unique  Green’s  function  for  the  air¬ 
bridge  structure  (Fig.  1).  This  was  accomplished  by  working  with  Maxwell’s 
equations  to  establish  a  tractable  equation  and  by  representation  of  our  source 
as  dirac  delta  functions.  We  then  applied  boundary  conditions  to  solve  for  the 
unknown  coefficients  and  formulated  a  solution.  We  now  have  expressions  that 
will  give  the  resulting  field  produced  by  a  point  source  directed  in  the  x  or  z 
direction  as  required  to  analyze  the  air-bridge  structure. 

2.3  Application  of  Method  of  Moments 

The  method  of  moments  is  a  numerical  technique  used  for  solving  functional 
equations  which  cannot  be  solved  in  closed  form.  By  reducing  the  functional 
relation  to  a  matrix  equation,  known  techniques  can  be  used  to  solve  the  resulting 
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I 


A 


matrix  equation.  This  method  is  computationally  intensive  but  with  the  advent 
of  faster  computers,  the  method  has  become  feasible. 

To  apply  the  method  of  moments  in  the  specific  case  of  an  air-bridge,  one  should 
follow  the  steps  outlined  below: 

1.  Use  the  integral  equation  (11)  derived  in  section  1  along  with  the  relations 
(13)  and  (5)  so  one  obtains  an  integral  equation  that  relates  the  current  to  the 
electric  and  magnetic  fields  respectively.  A  general  form  of  this  is 


(95) 


where  Lop  is  an  integral  operator  operating  along  with  the  derived  Green’s  function 
and  g  is  a  vector  function  of  either  the  electric  field  or  magnetic  field  H. 

2.  Represent  the  current  on  the  conducting  strip  aa  a  sum  of  coefficients  mul¬ 
tiplied  by  a  pre-determined  basis  function, 

Nia 

?.  =  £/, 7,  (96) 

qml 

In  equation  (96)  where  /,  represents  the  complex  coefficients,  Nia  represents  the 
number  of  sections  the  strip  is  to  be  divided  into  and  7,  represents  the  chosen 
basis  functions  which  represent  the  current  distribution. 

3.  Discretize  the  integral  equation  by  minimizing  the  resuting  error  function 
6E  on  the  surface  of  the  conducting  strips. 

In  applying  the  above  steps  to  the  problem  of  an  air-bridge  we  have  set  up  the 
integral  equation  for  the  electric  field  of  the  form 

T  =  ->u//i  f  j  J  (7  +  iw)  •  7?  •  IdV.  (97) 
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a.  One  dimensional  view  highlighting  current  directions 

Figure  6:  Currents  are  assigned  variables  by  direction. 

In  proceeding  to  the  second  step,  we  seperate  the  problem  into  five  different 
sections  as  pictured  below  (this  figure  is  taken  in  part  from  fig. I). 

For  J\x  and  we  will  model  the  current  as  a  sum  of  an  incident  current  A, 

a  reflected  current  B,  and  the  sum  of  incremental  currents  For  Jix,  J\,,  and 
Ju  only  a  sum  of  incremental  currents  is  required.  Implementing  this  convention 
results  in 

s\ * 

J„  =  (98) 

<J=l 

Nix 

Ju  —  ^  (^9) 

9*1 

,V3* 

h,  =  Ai,e>k*  +  Bi:e->kx  +  a*7**  (100> 

9*1 
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(101) 


Six 

Jlx  ~  y  (  IglxJ  1  za 

q=\ 


S2x  _ 

^2x  =  IglxJ 2xo- 

?=l 


(102) 


The  basis  functions  the  x-directed  currents  in  this  case  are  the  same  and  are  defined 
as  being 

(  ««"(*(*' _  c  T'  C  T 

~.)1  x  <  x  <  x 

nn(Ut)  Z9  -  1  ^  Z«+l 

<*v>) 


«W*\f)  - 


(-M - 1 - r  -*<«'<“ 

\*V)  /  ,  iX,\i  2  —  *  —  2 


(103) 


V 


(-(*)’)’ 


where  for  the  y-direction  we  have  used  a  maxwellian  distribution  function. 

The  first  step  in  discretizing  the  integral  equation  (97)  is  to  evaluate  the  error 
that  our  mathematical  representation  of  the  electric  field  will  produce.  This  is  done 
by  evaluating  the  electric  field  produced  by  one  section  of  current  on  the  conductor 
on  another  section  of  the  conductor.  Since  the  electric  field  on  any  part  of  the  strip 
must  be  zero,  the  value  of  our  integral  evaluates  the  error.  Then,  by  the  concept 
of  least  square  estimation,  when  one  takes  the  inner  product  of  the  basis  function 
along  with  the  error  and  sets  the  result  to  zero,  the  error  is  minimized 

Proceeding  to  do  this  for  our  problem  here  we  first  account  for  the  different 
components  of  the  fields  produced  by  a  current  of  given  orientation. 

Ex  —  E\gg  ■+■  Ei gg  +  Eifg  +  E\tt  +  fjr,  ( 104) 
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and 


(105) 


First,  we  shall  consider  the  x-directed  field  produced  by  a  x-directed  current  on 
the  dielectric 

£-  -  (t1’)  e~““dx'J-f  i-h)  r  _  (L)2iiJin  ^ dy 

p  tan(klzri(H  —  c)) 

ab(cos(kzH))[k[11  tan(kzH)  —  kttan(k[tl(H  —  c))] 


(*!  -  (t)1  cm  OfO sm  (rr)  + 

p  tan(*,ff)  (y)  (<»'  -  l)t<m(*{"(H  -  c)) 
a6co5(fc^)[fc'/rfan(fc,//')  -  kttan{k[lI{H  -  c))] 


-  c))  -  )1 


(?*)  (=)  (x >')  *' 


_ g  tan(k["{H  -  c)) _ 

oi(a»(*,ff))[*'"ian(t,W)  -  k,ton(k{"(H  -  c))| 


»*)  + 
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abcos(ktH)[k‘Iltan(k,H)  —  k,tan{klxn(H  —  c))] 


[kyltan(k*lI(H  —  c))  —  t*llkgtan(ktH)] 


(t)  (^r) Jin  (^r2)  cw(fc'z). 


tha-Y.rnmpnufl  pf  tU  .Urtrin  flalrl  nn  »^«>  ntl.^ 

For  a  x-directed  current,  a  z-directed  field  is  also  produced.  Proceeding  with 
the  same  procedure  as  above  one  obtains 


_  P  tan(kjn(H  -  c)) _ 

ab(cos{ktH))[kltIrtan(kxH)  -  kttanik^^H  -  c))] 


/—mr\  .  fntrx\  ( miry\ 

{—) sm  ItJ  Jm  {—)  + 

P  t<rn(k,H)  (y)  (4"  -  1  )tan(kj"(H  -  c)) 
abcos{kxH)[k[lltan{ktH)  —  kttan{kltll(H  -  c))] 

_ 1 _ 

[k[lItan(k[u(H  -  c))  -  t{.11  kttan(ktH)] 


( k 3  -  k])sin  sin  «n(fc,z)j  + 

L\  “*  (t*')  (£)  ^7^jTJ,n  (“5*# )  ^ 

' _ y  tan(k!,u{H  -  <Q) _ 

.a6(co3(lb1//))[Jb"^an(lbI^)  -  k,tan(k‘"(H  -  c))] 


/— wr\  /nirx\  /miry\  , 

{—) 3tn  {—) jw  rr ) 


■)+ 
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p  tan(k,H)  (“)  (c i"  -  l)tan(k{"(H  -  c)) 
abcos{k,H)[kilItan{k,H)  -  kxtan{kiIl{H  -  c))) 

1 

[Jfc'"tan(Jfc'"(tf  -  c))  -  elrllkxtan(kxH)] 


f  (-) 

J-*  \xwJ 


(mx  A  ,  » 


[*-(*) 

?  taw(fcf/z(g  -  c)) 


a6(co3(it,/f))[ik",<fln(ib,tf )  -  -  c))] 


t>  tan(k'H)  (y)  -  1  -  e)) 

a6cM(fc,//)[l;',,tan(fcs/f)  -  kttan(kltll(H  -  c))] 


1 

[klglltan(kltn{H  -  c))  - 
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( A: 1  —  k1)sin  sin  sin(kzz) 

For  a  z-directed  current,  one  must  consider  the  fields  above  and  below  the 
source  seprately  as  they  have  different  forms  in  those  two  respective  regions.  Also, 
as  before,  both  x-  and  z-components  of  the  electric  field  are  produced. 

First  we  consider  the  x-directed  field  produced  by  a  z-directed  current  above 
the  source  where 


(107) 


Eu 


iV,. 

= 

<1*1 


(v- 


ji’,(^'('.^'1>)(e;"t.coa(t.(r'  -  H))co,(k!"h)+ 


kilIsin(k,(z  -  H))sin(kilIh)dz)  + 


k[u  sin(kt(z  —  H))sin{k[tl h)dz)J 


4 ft _ 1 _  .  /  ntr  A 

ab  tlrnk\sin(ktH)co3(k[llh)  +  kltnkxcos{ktH)sm{kltIlh)Sm  V  a  *  ) 

(==)  CO,  (52)  *n  (23)  (-*.)»»<M)  (108) 

For  the  x-directed  field  below  the  source  produced  by  a  z-directed  current,  one 
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obtains 


'16r* 


=  f/„,  ( r  +  r* 

0=1  V*'*  sm{klx)  J o 


sm(klx) 


cos(kxz  )dz 


L\  (£)  — "n  (tt*)  ^ 


[«-(¥)T 


4^  eIrnkxsin(kx(z  -  H))cos(k[11  h)  4-  k^cosjk^z  —  //))sm(fcf;//t)  ^nrr 
kxeln3in(kxH)cos(klxtlh)  +  k^nsin(klnh)cos(kxH) 


■am 


( nir\  /nirx\  /mirtA 

It;  co3  (— ; 

For  the  z-directed  field  above  the  source  produced  by  a  z- directed  current,  one 
obtains 

£....  =  £  l,u  (  jf  -  H))co,(k[" k)+ 


(109) 


kjrisin(kx(z'  —  H))ain(klM11  h)dz  )  + 


-  H))ain(k["k)d,)) 
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4/i  1 

ab  elllk^sin{k,H)cos(k‘IIh)  +  k[u ktcos{k,H)sin(kIx11  h) 


( k 2  —  kj)ain  s*n  (~^)  cos(ktz) 


stn 


(¥*) 


(no) 


For  the  z-directed  field  below  the  source  produced  by  a  2-directed  current,  one 
obtains 

-  ,5  '"U  <in(Ifcj C0‘(k-Z  )iZ  +  l  - 7MM7) - ms,-k-z)iz 

ff  (±) - i - T»,,  (21y)  dy 

■'-f  \rwJ  ft  _  ^21*  \  b  y  J  y 


4/i  .  (nir  ,\ 

yn(  ~x) 


ejl 1  ktsin(kt(z  -  H))cos(k[n  h)  4-  kltncos(kt(z  —  H))sin(kltn  h) 
kMelr11  sin(ktH)co3{k,t11  h)  +  kltu 3in{klMn h)cos{k,H) 


(*  _  kD„n  (Sli)  Jin  (22) 


For  a  x-directed  field  in  the  region  above  a  x- directed  current  which  is  above 
the  dielectric,  one  obtains 
N u 

E 

*•1 


(111; 


'  larra 


£/„, ( /*  (six) 

ft  ,fc  V-,  «n(M,)  V  a  ) 


[*  I!S<*!£i£Zi]iCM  (Six')  ix  ) 

/o  s»n(e/,)  \  a  J  j 
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n  (-L)  - i -  3in  ( dy 

J-1  \irwJ  ^  b  > 

( k j  _  p  -  cot(k,H) ] 

V  V  a  )  )  abkt[cot(kMz')cos(ktz')  —  sm(/b,.z')] 


cos 


/ nirx 

V  a 


\  (mirx 

)  Stn  [— 


J  sin(ktz) 


p(l  -  tzn^k.H))  (M)  -  l)ian(*»'(ff  -  e» 

a6[^//<an(t,/f)  -  kttan{k[ll(H  -  c))][coi(fctz')coj(Jfc,z')  -  sin{ktz')\ 


k[n tan{k[n [H  -  c))  —  tilIkxtan(kxH) 


(Z)  k.cos  (==)  (=p)  (-k.s,n(k.z) 


(112) 


For  the  x-directed  field  located  below  the  x-directed  current  source  which  is 
above  the  dielectric,  we  obtain 

r  v*  r  (  f°  sin(k(z'  —  x,_t))  fnx  ,\  ,  ,  . 

£“'~  -  g '«•  [L,  ■  zi„(ki.)  ^{Tx)dx  + 


L 


!  sin(k(zt+x  - 
sin(kl 


(t1  )  dx) 


r  (i.) - ! - j,m  (==,')  dy 

r,_  (*)’]*  v  6  > 

f  ,  /n*\:\ _ f  [*"'  ~  k.tanik'J'm  -  c))tan(k,H) 

\  \  o  /  /  aWk,[cot(ibIz')coj(fc,z')  —  )][Arf 


')}[klIltan{k,H)  -  kttan{k‘!l i 
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fnirx\  .  i'miry\  .  . 

cos  —  J  stn  l — — J  stn(k,  z)  + 


_ —p _  fmrx\  frn‘xy\  ,.rr 

abk1[cot(ktz)cos(kxz')  —  Jtn(Ar,z')]C°3  \  a  )  5m  V  b  )  C°9  z  ^  + 

Ml  -  tan^t.g))  (*jf)  (t"'  -  l)tan(k‘"(H  -  cj) 
ab[k[ntan(kxH)  -  kxtan{k[u{H  -  c))j[cot(fczz')cos(fczz')  -  sin(fczz')] 

_ 1 _ 

k[utan{k[lI{H  -  c))  —  elrnkxtan(kxH) 

(t)  (^M^)  (-*.)«»(*.*) 

For  a  z-directed  field  above  the  current  produced  from  a  x-directed  current 
above  the  dielectric,  one  obtains  the  expression 
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sin  (—p-'j  ( k,)cos(k,z )  + 


P(  1  -  tan*(ktH))  (m)  (e1"  -  l)tan(ki“(H  -  c)) 
ab[k[lItan(kzH)  —  kztan(k[11  (H  —  c))][cof(fcz  z')cos(kzz')  —  ain(Jkzi')] 


kllitan(klIl(H  —  c))  —  tIrlIkztan(kzH) 


( fc 2  -  fc*)sin  sin  (~p)  c^(fc»2) 


(114) 


For  a  z-direct^d  field  below  the  current  produced  from  a  x-directed  current 
above  the  dielectric,  one  obtains  the  expression 


■  !'•*(/: 


1  sin(k(x  ~  fnir 

f  sin(klx)  003  V  a 


/nr  A  , 
cos  ^ — x  J  dx  + 


ri  sin(fc(g,+i  - 
3  sin(klx) 


-x'))  fnx  ,\  ,  ,\ 

y — «•  {tx  ) dx ) 


/*  (-?-) - i - r«n  (2=,')  dy 

J">  {rwl  [l -(*)’]*  U  } 


_ P  [kj!t  -  k,tan(kjIl{H  -  c))<an(fc,/f)] _ 

a6ife,(cot(fc,x')cos(ifciz')  -  sm(Jfc,z')][Jb£,;tan(fc,/f)  -  kztan{k[lI{H  -  c))] 


(~t)  (=r) J,n  (~p)  + 


L,  —  ,  <L  „  (-22)  3>n  (— )  sm  ™ 

l?f[co<(fc,z  )cos(ife,z  )  -  stn(kzz  ))\  a  J  \  a  J  \  b  J 
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_ p{  1  -  tan*(ktH))  («)  (ej11  -  l)tan(kjn(H  -  c)) _ 

a6[A:///tan(fc,/f)  -  kztan{k[n{H  -  c))][cot(fc,z')cos(fc,z')  -  sin(^z')] 

_ 1 _ 

klxIltan(klzlI(H  —  c))  —  ef.11  kxtan(ktH) 

( k 3  —  k])sin  sin  cos(kxz) 

2.4  Matrix  Equation 

The  resulting  matrix  equation  is  formed  and  upon  inversion  the  unknowns  can 
be  obtained  for  a  given  exciation. 

3  Scattering  Parameters 

Using  the  derived  current  distribution  on  the  conductors,  one  can  apply  an  ideal 
transmission  line  model  to  determine  the  scattering  parameters. 

4  Summary 

From  this  general  analysis  of  the  Green’s  function  and  the  ultimate  detei  na¬ 
tion  of  the  scattering  parameters  of  an  air-bridge,  the  formulation  used  here  could 
be  applied  to  a  variety  of  structures  whose  geometry  requires  a  three-dimentional 
analysis. 

Other  planned  work  includes  the  application  of  this  work  to  other  structures 
and  the  use  of  air-bridges  as  a  curcuit  element  in  the  construction  of  other  passive 
microwave  circuits. 
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